New evidence for super-roughening in crystalline surfaces with disordered substrate

We study the behavior of the Binder cumulant related to long distance correlation functions of the discrete Gaussian model of disordered substrate crystalline surfaces. We exhibit numerical evidence that the non-Gaussian behavior in the low-$T$ region persists on large length scales, in agreement with the broken phase being super-rough.Comment: 10 pages and 4 figures, available at http://chimera.roma1.infn.it/index_papers_complex.html . We have extended the RG discussion and minor changes in the tex

Ground state properties of solid-on-solid models with disordered substrates

We study the glassy super-rough phase of a class of solid-on-solid models with a disordered substrate in the limit of vanishing temperature by means of exact ground states, which we determine with a newly developed minimum cost flow algorithm. Results for the height-height correlation function are compared with analytical and numerical predictions. The domain wall energy of a boundary induced step grows logarithmically with system size, indicating the marginal stability of the ground state, and the fractal dimension of the step is estimated. The sensibility of the ground state with respect to infinitesimal variations of the quenched disorder is analyzed.Comment: 4 pages RevTeX, 3 eps-figures include

Using network-flow techniques to solve an optimization problem from surface-physics

The solid-on-solid model provides a commonly used framework for the description of surfaces. In the last years it has been extended in order to investigate the effect of defects in the bulk on the roughness of the surface. The determination of the ground state of this model leads to a combinatorial problem, which is reduced to an uncapacitated, convex minimum-circulation problem. We will show that the successive shortest path algorithm solves the problem in polynomial time.Comment: 8 Pages LaTeX, using Elsevier preprint style (macros included

Dislocations in the ground state of the solid-on-solid model on a disordered substrate

We investigate the effects of topological defects (dislocations) to the ground state of the solid-on-solid (SOS) model on a simple cubic disordered substrate utilizing the min-cost-flow algorithm from combinatorial optimization. The dislocations are found to destabilize and destroy the elastic phase, particularly when the defects are placed only in partially optimized positions. For multi defect pairs their density decreases exponentially with the vortex core energy. Their mean distance has a maximum depending on the vortex core energy and system size, which gives a fractal dimension of $1.27 \pm 0.02$. The maximal mean distances correspond to special vortex core energies for which the scaling behavior of the density of dislocations change from a pure exponential decay to a stretched one. Furthermore, an extra introduced vortex pair is screened due to the disorder-induced defects and its energy is linear in the vortex core energy.Comment: 6 pages RevTeX, eps figures include

Super-roughening as a disorder-dominated flat phase

We study the phenomenon of super-roughening found on surfaces growing on disordered substrates. We consider a one-dimensional version of the problem for which the pure, ordered model exhibits a roughening phase transition. Extensive numerical simulations combined with analytical approximations indicate that super-roughening is a regime of asymptotically flat surfaces with non-trivial, rough short-scale features arising from the competition between surface tension and disorder. Based on this evidence and on previous simulations of the two-dimensional Random sine-Gordon model [Sanchez et al., Phys. Rev. E 62, 3219 (2000)], we argue that this scenario is general and explains equally well the hitherto poorly understood two-dimensional case.Comment: 7 pages, 4 figures. Accepted for publication in Europhysics Letter

Numerical study of the strongly screened vortex glass model in an external field

The vortex glass model for a disordered high-T_c superconductor in an external magnetic field is studied in the strong screening limit. With exact ground state (i.e. T=0) calculations we show that 1) the ground state of the vortex configuration varies drastically with infinitesimal variations of the strength of the external field, 2) the minimum energy of global excitation loops of length scale L do not depend on the strength of the external field, however 3) the excitation loops themself depend sensibly on the field. From 2) we infer the absence of a true superconducting state at any finite temperature independent of the external field.Comment: 6 pages RevTeX, 5 eps-figures include

Ground state properties of fluxlines in a disordered environment

A new numerical method to calculate exact ground states of multi-fluxline systems with quenched disorder is presented, which is based on the minimum cost flow algorithm from combinatorial optimization. We discuss several models that can be studied with this method including their specific implementations, physically relevant observables and results: 1) the N-line model with N fluxlines (or directed polymers) in a d-dimensional environment with point and/or columnar disorder and hard or soft core repulsion; 2) the vortex glass model for a disordered superconductor in the strong screening limit and 3) the Sine-Gordon model with random pase shifts in the strong coupling limit.Comment: 4 pages RevTeX, 3 eps-figures include

Application of a minimum cost flow algorithm to the three-dimensional gauge glass model with screening

We study the three-dimensional gauge glass model in the limit of strong screening by using a minimum cost flow algorithm, enabling us to obtain EXACT ground states for systems of linear size L<=48. By calculating the domain-wall energy, we obtain the stiffness exponent theta = -0.95+/-0.03, indicating the absence of a finite temperature phase transition, and the thermal exponent nu=1.05+/-0.03. We discuss the sensitivity of the ground state with respect to small perturbations of the disorder and determine the overlap length, which is characterized by the chaos exponent zeta=3.9+/-0.2, implying strong chaos.Comment: 4 pages RevTeX, 2 eps-figures include

Shunting operations at flat yards : retrieving freight railcars from storage tracks

In this paper, we study the railcar retrieval problem (RRT) where specified numbers of certain types of railcars have to be withdrawn from the storage tracks of a flat yard. This task arises in the daily operations of workshop yards for railcar maintenance. The objective is to minimize the total cost of shunting via methods such as minimizing the usage of shunting engines. We describe the RRT formally, present a mixed-integer program formulation, and prove the general case to be NP-hard. For some special cases, exact algorithms with polynomial runtimes are proposed. We also analyze several intuitive heuristic solution approaches motivated by observed real-world planning routines. We evaluate their average performances in simulations with different scenarios and provide their worst-case performance guarantee. We show that although the analyzed heuristics result in much better solutions than the naive planning approach, they are still on average 30%-50% from the optimal objective value and may result in up to 14 times higher costs in the worst case. Therefore, we conclude that optimization should be implemented in practice in order to save valuable resources. Furthermore, we analyze the impacts of yard layout and the widespread organizational routine of presorting on the railcar retrieval cost

The critical exponents of the two-dimensional Ising spin glass revisited: Exact Ground State Calculations and Monte Carlo Simulations

The critical exponents for $T\to0$ of the two-dimensional Ising spin glass model with Gaussian couplings are determined with the help of exact ground states for system sizes up to $L=50$ and by a Monte Carlo study of a pseudo-ferromagnetic order parameter. We obtain: for the stiffness exponent $y(=\theta)=-0.281\pm0.002$, for the magnetic exponent $\delta=1.48 \pm 0.01$ and for the chaos exponent $\zeta=1.05\pm0.05$. From Monte Carlo simulations we get the thermal exponent $\nu=3.6\pm0.2$. The scaling prediction $y=-1/\nu$ is fulfilled within the error bars, whereas there is a disagreement with the relation $y=1-\delta$.Comment: 8 pages RevTeX, 7 eps-figures include